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We study the Bethe-Salpeter equations for spin zero diquark composites in the color superconducting phase 
of Nf = 2 cold dense QCD. The explicit form of the spectrum of the diquarks, containing an infinite tower of 
narrow (at high density) resonances, is derived. It is argued that there are five pseudo-Nambu-Goldstone bosons 
(pseudoscalars) that remain almost massless at large chemical potential. These five pseudoscalars should play an 
important role in the infrared dynamics of Nf = 2 dense QCD. 



In his studies Dmitrij Vasilievich Volkov was 
always led by the beauty within the problems he 
considered. One of his passions was the theory 
of spontaneous symmetry breaking, in particular, 
the dynamics of Nambu-Goldstone (NG) parti- 
cles (both bosons and fermions), to which he con- 
tributed a great deal to our present understanding 
|IJ|^]. Therefore we think it is most appropriate 
to report in this volume a recent investigation of 
the dynamics of diquark composites (in particu- 
lar, diquark NG bosons) in cold dense QCD with 
two fcrmion flavors. 

Only a few years ago, not much was known 
about the properties of different phases in dense 
quark matter (see, however, Refs. H,|[). The situ- 
ation drastically changed after the ground break- 
ing estimates of the color superconducting order 
parameter were obtained in Refs. HBl. Within 
the framework of a phenomenological model, it 
was shown that the order parameter could be as 
large as 100 MeV. Afterwards, the same estimates 
were also obtained within the microscopic theory, 
quantum chromodynamics fi-[lj]. The further 
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progress in the field was mostly motivated by the 
hope that the color superconducting phase could 
be produced either in heavy ion experiment, or in 
the interior of neutron (or rather quark) stars. 

Despite many advances |lj||l^] in study of the 
color superconducting phase of dense quark mat- 
ter, the detailed spectrum of the diquark bound 
states (mesons) is still poorly known. In fact, 
most of the studies deal with the NG bosons of 
the three flavor QCD. At best, the indirect meth- 
ods of Refs. rp~5| — [3~sfl could probe the properties of 
the pseudo-NG bosons. It was argued in Ref. fll9f , 
however, that, because of long-range interactions 
mediated by the gluons of the magnetic type , 
the presence of an infinite tower of massive di- 
quark states could be the key signature of the 
color superconducting phase of dense quark mat- 
ter. 

In this paper, we consider the problem of spin 
zero bound states in the two flavor color supercon- 
ductor using the Bethe-Salpeter (BS) equations. 
We find that the spectrum contains five (nearly) 
massless states and an infinite tower of massive 
singlets with respect to the unbroken SU(2) C sub- 
group. Furthermore, in the hard dense loop im- 
proved ladder approximation, the following mass 
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formula is derived for the singlets: 

where k is a constant of order 1 (we find that 
k ~ 0.27), |A| is the dynamical Majorana mass 
of quarks in the color superconducting phase, and 
a s is the value of the running coupling constant 
related to the scale of the chemical potential \i. 

At large chemical potential, we also notice an 
approximate degeneracy between scalar and pseu- 
doscalar channels. As a result of this parity dou- 
bling, the massive diquark states come in pairs. 
In addition, there also exist five massless scalars 
and five (nearly) massless pseudoscalars [a dou- 
blet, an antidoublet and a singlet under SU(2) C ). 
While the scalars are removed from the spec- 
trum of physical particles by the Higgs mecha- 
nism, the pseudoscalars remain in the spectrum, 
and they are the relevant degrees of freedom of 
the infrared dynamics. At high density, the mas- 
sive and (nearly) massless states are narrow res- 
onances. 

In the case of two flavor dense QCD, the orig- 
inal gauge symmetry SU (3) c breaks down to 
the SU(2) C by Higgs mechanism. The flavor 
SU(2) l x SU(2)r group remains intact at the vac- 
uum. The appropriate order parameter is an an- 
titriplet in color and a singlet in flavor. Without 
loss of generality, we assume that the order pa- 
rameter points in the third direction of the color 
space. In order to have a convenient description 
of the bound states at the true vacuum, we intro- 
duce the following Majorana spinors, 

K=ri+^ ij &% a =1,2, (2) 
*a = C-£3a 6 e y W> C )?, a =1,2, (3) 
made of the Weyl spinors of the first two colors, 

1>i = V + (*D)i, ^ c ))=V-^ c D ) b v (4) 

<t>i=V-(*D)l {4> c ) h j =V + (9^) b j . (5) 

Here i,j = 1,2 are flavor indices, V± — (l±7 5 )/2 
are the left- and right-handed projectors, $ d is 
the Dirac spinor, and ^£ = is its charge 

conjugate. Regarding the quark of the third color, 
we use the Weyl spinors, ip l and <fi l , for left and 



right components, respectively (notice that the 
color index is omitted). 

The BS wave functions of the bound diquark 
states in the channels of interest are given by 

(2tt) 4 6 4 (p + - P --P)x^(p,P) = 

= (o\T¥ a ( P+ )$ i (-p-)\P;b) L , (6) 

(2ny6H P+ -P--P)\ h a ) (p,P) = 

= (0|T^(p + )*J(-p_)|P;5) £ , (7) 
(2w) 4 5 i ( P+ -p--P)r,(p,P) = 

= (0\T¥ a (p + )^i-P-)\P)L, (8) 
(2*)*S i (p+-p--P)*(p,P) = 

= (0\T^( P+ )M-p_)\P) L , (9) 

where p = (p + + p_)/2 and the quantities on the 
right hand side of these equations are defined as 
the Fourier transforms of the corresponding BS 
wave functions in the coordinate space. There are 
also the BS wave functions constructed out of the 
right handed fields Q l a and cj> % . One might notice 
that there is another diquark channel, a triplet 
under SU (2) c , that we do not consider here. The 
reason is that the repulsion dominates in such a 
channel, and no bound states are expected [the 
triplet comes from the SU(3) C sextet]. 

In order to derive the BS equations, we use 
the method developed in Rcf. ^(J for the case of 
zero chemical potential. To this end, we need to 
know the quark propagators and the quark-gluon 
interactions. 

By introducing the multicomponent spinor that 
combines the Majorana spinors of the first two 
colors and the Weyl spinors of the third color, 

ipf^j , we find that the inverse propaga- 
tor takes the following block-diagonal form: 

G- 1 = diag , s; 1 ^- , s- 1 8/) , (10) 

where, upon neglecting the wave functions renor- 
malization of quarks |q-[T^1 , 

S- 1 = -*(V+M7°7 5 + A p P-+A p P + ) , (II) 
8- 1 = -i(j + vrf > )V+, (12) 

a; 1 = -i{j-i*f)r-. (13) 

Here the notation, A p = A+A+ + , A p = 



3 



7°A p 7 C \ and A^ = (1 ±5 • are the same 

as in Ref . || . 

The bare vertex, 7^, is also a 3 x 3 matrix, 

732% 733V 



where 
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721° j 



7lV 
723^ 



(14) 



with 



(71!): 



T^ 6 - 25 A T* b V- 



I22 
7 3 A 3 



_jiA3 



V-, 



(7^2) 
(7&) 
(7 21 ) 
(731) 

(73i) 



T A3 V+, 

-e 3ac T 3 Ac V- 

T 3 Ab V+, 

0, 
0. 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 

(23) 



where T A are the SU (3) c generators in the funda- 
mental representation. By making use of this ver- 
tex and the propagator in Eq. (|l0|), it is straight- 
forward to derive the BS equations in the (hard 
dense loop improved) ladder approximation. The 
details of the derivation, as well as the explicit 
form of equations are given elsewhere pj] . Here 
we just note that the most transparent form of 
the equations appears for the amputated BS wave 
functions, defined by 



xip,P) = 


S -1 &H 
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In order to get a feeling of the problem at hand, 
let us briefly discuss the analysis of the BS equa- 
tion for the %-doublet. In general, the BS wave 
function contains eight different Dirac structures 



22 j . It is of great advantage to notice that only 
four of them survive in the center of mass frame, 
P=(M b ,0), 



A. a 



{b) (P,0)=SaX(p), 



(28) 



xip) 



Xi A+ + (po 



S + -7r)X2 7 A+ 



-xtK + (po + 4 + ^)xh°K 



(29) 



with €p — \p\±n [the factors (po±€p + Mb/2) are 
introduced here for convenience] . This is the most 
general structure that is allowed by the space- 
time symmetries of the model. 

Now, in the particular case of the NG bosons, 
Mf, = 0, we will show that the BS wave function 
is fixed by the Ward identities. Indeed, let us 
consider the following non-amputated vertex: 

r#„(*.tf) = (0\Tj A (0)¥ a (x)^(y)\0), (30) 

where, for our purposes, it is sufficient to consider 
A = 4, . . . , 8 (that correspond to the five broken 
generators). In the (hard dense loop improved) 
ladder approximation, the vertex satisfies the fol- 
lowing Ward identity [Ml : 



P^aLik- 



P,k) 



iT A3 5)[ Sk -S k+P }T-. (31) 



As in the case of the BS wave functions, it is 
more convenient to deal with the corresponding 
amputated quantity, 



^fJk + P,k) 



0-1 w 

J k+P 3 - aj,ti\ 



(32) 



(k + P,k)s^. 

This latter satisfies the following identity: 

P<t(k+P k) = iT A H) [S£ P - s?\ 7> + .(33) 

By making use of the explicit form of the quark 
propagators in Eqs. (|ll|) and (|l^), we could check 
that the right hand side of Eq. (|3^) is non-zero 
in the limit P — > 0. This is possible only if the 
vertex on the left hand side develops a pole as 
P — > 0. After a simple calculation, we obtain 



KL(k + Pk) 



p py « 3 



(34) 



where, we introduced P M = (Po, c^P) with c x be- 
ing the velocity of the NG boson in the x-doublet 
channel. 

By making use of the definition in Eqs. (||) and 
(p4|), it is also not difficult to show that the pole 
contribution to the vertex function (53) is directly 
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related to the BS wave function. By omitting the 
details, 

X £%0) = S^ x (p,0) = 6^^V + , (35) 

where is the decay constant of the corre- 

sponding doublet whose formal definition is given 
by 



(0\^2T^(0)\P,b) L = iSlP^. 



(36) 



A=A 



By comparing the Dirac structures in Eqs. ( J29| ) 
and p5]), we see that no components of the \i 
type appear in Eq. (|3^) which follows from the 
Ward identities. It was rewarding to establish 
that, in this approximation, the structure of the 
BS wave function required by the Ward identity 
is indeed a solution to the BS equation for the 
X-doublet. A similar situation takes place for the 
77-singlet p^| . 

Now let us discuss the fate of the massless 
states that we obtain. Altogether, there are five 
scalars and five pseudoscalars (a doublet, an an- 
tidoublet and a singlet). Because of the Higgs 
mechanism, the scalars are removed from the 
spectrum. Nevertheless, these scalar bound states 
exist in the theory as "ghosts" |^4j , and one can- 
not get rid of them completely, unless a unitary 
gauge is found. In fact, these ghosts play a very 
important role in getting rid of unphysical poles 
from the on-shell scattering amplitudes p4| ]. 

As for the pseudoscalars, they remain in the 
spectrum as pseudo-NG bosons. In the (hard 
dense loop improved) ladder approximation, they 
look like NG bosons because the left and right 
sectors of quarks decouple. One could think of 
this as an effective enlargement of the original 
color symmetry from SU (3) c to an approximate 
SU(3)cx x SU(3) C} r. Then, since the approxi- 
mate symmetry of the ground state is SU (2) c .l x 
SU(2) Ct R, five scalar NG bosons (which are re- 
moved by the Higgs mechanism) and five pseu- 
doscalar NG bosons (which remain in the spec- 
trum) should appear. Of course, in the full 
theory, the pseudoscalars are only pseudo-NG 
bosons. Indeed, they should get non-zero masses 
due to higher orders corrections that are beyond 



the improved ladder approximation fl25j . Since 
the theory is weakly coupled at large chemical 
potential, it is natural to expect that the masses 
of the pseudo-NG bosons are small compared to 
the value of the dynamical quark mass. 

We conclude our discussion of the massless di- 
quarks by emphasizing that the low-energy dy- 
namics of the two flavor QCD is dominated by 
massless quarks of the third color (which might 
eventually get a small mass too if another (non- 
scalar) condensate is generated |i|26|) and by the 
five pseudoscalars that remain almost massless in 
the dense quark matter. Of course, the gluons 
(glueballs) of the unbroken SU(2) C may also be 
of some relevance but we do not study this ques- 
tion here. 

Now, let us consider massive diquarks. The 
structure of the BS equations becomes even more 
complicated in this case. In addition, one does 
not have a rigorous argument to neglect the com- 
ponent functions like xt m Eq. (|29|). In spite of 
this, we argue that all the approximations made 
before might still be reliable. Indeed, from the 
experience of solving the gap equation (which 
coincides with the BS equation for the mass- 
less states), we know that the most important 
region of momenta in the integral equation is 
A I <C p <C M- In this region, the kernel of the 
BS equations for massive states, A/&< |A|, is al- 
most the same. The deviations appear only in the 
infrared region where p< |A|. 

Therefore, in our analysis of the BS equations 
for massive states, we closely follow the approxi- 
mation used for the massless diquarks. By assum- 
ing that the component functions depend only on 
the time component of the momentum (compare 
with the analysis of the gap equation in Refs. 
|l4|]), we arrive at the following equation for the 
BS wave function of the massive singlets: 



?7i 0) = 



4tt 



dqK^(q)^(q)\n- 



A 



where A = (^f/ 2 ii/aj 2 
K^)( q ) = 



Al 2 



(37) 

and the kernel reads 
(38) 



g2+|A|*-(M„/2y" 

[Vi (p) is a scalar function that appears in the 
decomposition of the BS wave function 77 over 
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the Dirac matrices (compare with Eq. (f29|)]. At 
this point it is appropriate to emphasize that the 
Meissner effect plays an important role in the 
analysis of the massive bound states. Indeed, 
our analysis shows that these massive states are 
quasiclassical in nature, i.e., their binding energy 
is small compared to the value of the gap [see 
Eq. ([!])]. As a result, only the long range in- 
teraction mediated by the unscreened gluons of 
the unbroken SU(2) C is strong enough to produce 
these diquark states. We took this into account 
in Eq. ([37]). For completeness, we mention that 
the (nearly) massless (pseudo-) NG bosons are 
tightly bound, and the Meissner effect is not so 
important for their binding dynamics. 

By approximating the kernel ( |38| ) in each 
of the following three regions: < q < 

^\A\*-{M n /2)\ ^/|AP- (AV2) 2 < q < \A\ 
and A | < q < A, we could solve the BS equation 
( j37| ) analytically. Then, by matching the loga- 
rithmic derivatives of the separate solutions, we 
obtain the spectrum of the massive diquarks. By 
omitting the details, it is presented in Eq. ([j]). 

We would like to emphasize that for large // 
the hard dense loop improved ladder approxima- 
tion is reliable for the description of those bound 
states. The point is that a) the region of mo- 
menta primarly responsible for the formation of 
these composites is Eund%j Q *C /i, where the 
binding energy E bind ~ a 2 s A, and b) E bin d — ► oo 
as n — > oo §-0. Therefore the vacuum effects 
are higher order ones in a s in that region. Be- 
cause of that, the hard dense loop improved lad- 
der approximation, in which the contribution of 
the vacuum effects to the running of the coupling 
constant is neglected and only the running due 
to the polarization effects provided by the quark 
matter (non-zero fi) is taken into account, is jus- 
tifiable for large fi. 

Now, let us consider the case of massive di- 
quarks in the doublet channel. As is easy to 
check, the binding interaction in this channel is 
exclusively due to the five gluons affected by the 
Meissner effect. The approximate BS equation 
looks similar to Eq. (|3^), but with a different 
kernel and \q — p\ in the logarithm replaced by 
\q — p\ + cA where c = 0(1) is a constant 



At high density when the coupling constant is 
weak, this equation does not allow a non-trivial 
solution for M ^ 0. From the physical view- 
point, this indicates that the heavy gluons, with 
M g i ~ (a s /z 2 A) 1 / 3 3> A, cannot provide a suf- 
ficiently strong attraction to form massive radial 
excitations of the NG and pseudo-NG bosons. 

At the end, let us note that the massive diquark 
states may truly be just resonances in the full the- 
ory, since they could decay into the pseudo-NG 
bosons and/or gluons (glueballs) of the unbroken 
SU(2) C . At high density, however, both the run- 
ning coupling a s (p) and the effective Yukawa cou- 
pling g Y = \A\/F~\A\/n are small, 

and, therefore, these massive resonances are nar- 
row. 

In conclusion, in this paper we studied the 
problem of diquark bound states in the color 
superconducting phase of Nf — 2 dense QCD. 
While the scalar NG bosons are ghosts in the 
theory, the pseudoscalar pseudo-NG bosons are 
physical particles that should play an important 
role in the infrared. We also obtained the spec- 
trum of the massive narrow diquark resonances, 
whose existence would be a clear signature of the 
unscreened long range forces in dense QCD. 
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